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Possible Verification of Tilted Anisotropic Dirac Cone in a;-(BEDT-TTF)2 I3 Using 

Interlayer Magnetoresistance 
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It is proposed that the presence of a tihed and anisotropic Dirac cone can be verified 
using the interlayer magnetoresistance in the layered Dirac fermion system, which is realized 
in quasi- two-dimensional organic compound q;-(BEDT-TTF)2 I3. Theoretical formula for 
the interlayer magnetoresistance is derived using the analytic Landau level wave functions 
and assuming local tunneling of electrons. It is shown that the resistivity as a function of 
the azimuthal angle of the magnetic field takes the maximum in the direction of the tilt if 
anisotropy of the Fermi velocity of the Dirac cone is small. The procedure is described to 
determine the parameters of the tilt and anisotropy. 
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There is a group of condensed matter systems in which low-lying properties of the conduction electrons 
are described by a relativistic Dirac equation with the velocity of light replaced by the Fermi velocity. 
Remarkable physical phenomena of such Dirac fermions are clearly demonstrated in graphene, a monolayer 
of graphite.^' ^ In the observed integer quantum Hall effect the plateaus in the Hall conductivity are 
characterized by not integers, N ^ but N + 1/2. The origin of the shift of 1/2 is the presence of the zero 
energy Landau level for the Dirac fermions. 

In contrast to graphene, which is a purely two-dimensional system, the organic conductor a-(BEDT- 
TTF)2 Is^ is the first bulk material where the massless Dirac fermion-like spectrum is realized. Experi- 
mentally it was reported that this quasi-two-dimensional compound was a narrow gap system under high 
pressure. Using the tight-binding model with the transfer integrals obtained by X-ray diffraction ex- 
periment,^ Kobayashi et al. suggested that the system is zero-gap and the energy dispersion is linear 
around the zero-gap point. They showed that the electronic band structure is described by a tilted and 
anisotropic Dirac cone. The first principle calculations supported this Dirac cone structure. ^^'^^ Although 
observed negative interlayer magnetoresistance^^ is an evidence of the presence of a Dirac cone, as theo- 
retically explained by Osada^^ in terms of the zero energy Landau level wave function, there is no direct 
experimental verification that the cone is tilted and anisotropic. 

In this Letter, we propose that the Dirac cone structure is verified by analyzing interlayer magnetore- 
sistance through its dependence on the applied magnetic field direction. Using the analytic form of the 
zero-energy Landau level wave function with tilt and anisotropy of the Dirac cone,^^ we derive a formula 
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Fig. 1. (a)The definition of the crystal axes denoted by a and b and the wave vector axes and ky. The angle 
between the a-axis and the fca;-axis is represented by 0o. The arrow indicates the tilt direction. The angle between 
the fc^-axis and the tilt direction is defined as 7. By tilting the cone the intersection becomes an ellipse whose 
long axis is in the tilt direction as schematically shown in the figure. (b)The definition of the tilt angle, 9t. The 
horizontal axis is the tilt direction in the kx-ky plane. The vertical axis is the energy axis. The declined solid 
lines represent the dispersion of the Dirac cone along the tilt direction. The dotted line represents the center 
line of the two dispersions. The angle between the energy axis and the center line is 9t- 



for the interlayer magnetoresistance. Our formula is unique for the analysis of the tilted and anisotropic 
Dirac cone in a-(BEDT-TTF)2l3. It should be noted that in a-(BEDT-TTF)2l3 Shubnikov-de Haas os- 
cillations, which is a standard experiment to study an electronic structure, have never been observed and 
angle resolved photoemission spectroscopy, which directly confirmed the linear dispersion in graphene,^^ 
is not applicable to organic compounds. 

Before beginning the analysis of the system, we introduce parameters defining the tilt and anisotropy 
of the Dirac cone. We represent the crystal axes in the plane as a and b. We assume that the intersections 
of an anisotropic Dirac cone are elliptic. The principal axes are not necessarily parallel to the crystal axes. 
So the kx- and ky-axes in the wave vector space are taken so that those axes are parallel to the principal 
axes. The angle between the kx-axis and the a-axis is denoted by (pQ as shown in Fig. 1(a). As for the tilt 
direction, we introduce 7 to denote the angle between the k^-axis and the tilt direction. The parameter 
9t describes the tilt angle as defined in Fig. 1(b). Thus, the parameters are 7, 9t, and the Fermi velocity 
anisotropy (denoted by a below). Within our formulation, (j)Q is taken as a given parameter. 

Now we start with the Hamiltonian describing a single layer of an anisotropic and tilted Dirac cone 
system following Ref.l7, 

H=Y,H{kx,ky), (1) 

where 

H {kx, ky) = h( ''o^^ + ''o^^ ''f ^ " '""f' ] . (2) 

\ Vxkx -\~ iVyky Vq kx -\- VQky J 



This is a sufficiently general form. Anisotropy in the Fermi velocity is parameterized by a = Vx/vy. In 
the absence of the magnetic field the energy dispersion is given hy = h ( v^kx + v^ky it Jvxk"^ + f^/c^ j . 



2/8 



J. Phys. Soc. Jpn. Letter 
The angle 7 defined in Fig. 1(a) is given by 



7 = cot-^». (3) 



To compute the interlayer magnetoresistance, we need the Landau level wave functions. We represent 
the magnetic field as {Bx,By,Bz) = B {cos 9 cos (f), cos 6 sin (j), sin 9). Here cf) is the angle in the plane with 
respect to the positive kx-sxis and 9 is the angle between the magnetic field direction and the plane. 
We choose the gauge so that the vector potential is given by = ByZ + A^f^ Ay = —B^z + A^^^\ and 

(z) (z) 

A^ = where B^ = d^A];' - dyA'^^' . The presence of the inplane magnetic field is taken into account by a 
gauge transformation because they depend only on z. The zero-energy Landau level wave function is given 
in Ref.l5. Here we derive it in a different way which can be applicable to compute other Landau level 
wave functions. The derivation procedure consists of three steps. After the Pierls substitution {kx,ky) 
{kx, Ky) = {kx + eA^"^ /ch, ky + eA!y^ /ch), first we rescale ky so that the Fermi velocity in this direction is 
Vx- At this transformation, Bz is multiplied by Vy/vx- Secondly we rotate the system by the angle 7 in the 
plane. The transformed Hamiltonian is, 

H {kx, Ky) = hvxU^ {-riKxao + KxCTi + Kya2) U, (4) 

where Kx and Ky are redefined as the rotated variables, aj {j = 1, 2) are the Pauli matrices and 

(5) 

uJ^^i^^ ' ]. (6) 

V exp(-|7) J 

Now the direction of tilt is in the Kj;-axis. The angle of tilt is defined by 9t = """2 ~ , where tan 9± = l±rj. 
In the Schrodinger equation with the Hamiltonian UHU\ we subtract the term associated with the tilt 
from the both side of the equation. Applying the operator of the Dirac cone and after some algebra we 
obtain the operator form of an anisotropic harmonic oscillator. Diagonalizing the operator, we find the 
energies of the Landau levels. 




e„ = sgn(n) {hvx/lz) yj2\^\n\ (7) 



with n = 0, ±1, ±2, the magnetic length = ych/eBz and 



x\2 / y\2 



The Landau level wave functions are obtained by taking the Landau gauge for {A^c^ , Ay^"^ ) , 

/1,.|(CI"') . (9) 
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with sgn(O) = 1. Here k is the wave vector for the plane wave component, where the plane wave part is 
implicitly included, and 



^ (xsin7 + ycos7) + —k 



77A/2|n|sgn(n), (10) 



(-1)" „ ,,,_.J 1,2 



with Hn{C) the Hermite functions. 

The interlayer conductivity is calculated using the Kubo formula. As for the impurity scattering, we 
assume that the scattering leads to Lorentzian-shape density of states with half value width of a constant 
r. Although a self-consistent Born approximation shows that T is magnetic field dependent, here we 
ignore field dependence of F for simplicity. We focus on the contribution from the zero energy Landau 
level at zero temperature. We shall comment on the effect of the other Landau levels later. The matrix 
elements of the current operator is calculated similarly to the non-tilted case^"^ assuming local tunneling 
of electrons between two neighboring layers. The interlayer magnetoresistance is given by 



,{0) 



A 



P. 

Bq + B sin 6 exp 

where 



-hit) iS|^('^>«>7,A) 



(12) 



/((/), a, 7, A) = A ( a sin ^ cos 7 C0S(?!>sin7 

V a 

1 / 1 V 

Q;sin(^sin7-| — cose/) cos 7 . (13) 
A \ a ) 

Here Bq is a parameter of the theory and A = ^ ( '^n^ ) (^) taken to be a constant with 
the number of layers, tc the hopping parameter between neighboring layers, and Cc the lattice constant 
perpendicular to the plane. Note that the function /((^,q;,7, A) has the period tt with respect to <^ as 
confirmed analytically. Note also that there is another Dirac cone with 77 — —77 by symmetry.^^ But that 
Dirac cone has the same contribution because the expression is invariant under this sign change of i]. 

Let us move on to the parameter dependence of the interlayer magnetoresistance. Figure 2 shows the 
angular ^ dependence of ^p^} = —[p^zz — Po\/po for various values of A with a = l{vx = Vy). Here po is 
the interlayer resistance in the absence of the magnetic field. On the other hand. Fig. 3 shows the angular 
<p dependence for various values of a in the absence of the tilt. Anisotropy of the Fermi velocity also leads 
to an angular dependence of the interlayer magnetoresistance. As for the 9 dependence, we found that the 
ratio of the maximum to the minimum in Figs. 2 and 3 decreases with increasing 0. 

The fact that anisotropy in eq.(12) in the x — y plane arises from the tilt is understood as follows. 
In the presence of tilt, the intersection of the Dirac cone is deformed. In case of 77 = and Vx = Vy, the 
intersection is circle. But if 77 7^ 0, the intersection becomes an ellipse with the origin at {k^\ky^^) = 
{i]{en/hvx)/{l — '7^),0) and the ratio of the principal axes being A. The form of the Landau level wave 
functions is deformed according to this change of the intersection. We find that the Landau level wave 
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Fig. 2. Dependence of the interlayer magnetoresistance on (j) for various tilt parameters X a,t B = 6T, Bq = 0.5T 
and 6 = 20 degrees with a = l{vx = Vy). 
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Fig. 3. Dependence of the interlayer magnetoresistance on (p in the absence of the tilt (A = 1) for various a values. 
Other parameters are the same as in Fig. 2. 

function is more localized in the tilt direction than the perpendicular direction to the tilt direction. Since 
the matrix element of the interlayer current operator has a Fourier transformation like form with respect 
to the in-plane magnetic field, the current operator matrix element takes the minimum if the in-plane 
magnetic field is in the tilt direction. Therefore, the interlayer resistivity takes the maximum if the in- 
plane magnetic field is in the tilt direction. 

Kobayashi et al. estimated the parameters of the anisotropic tilted Dirac cone at the uniaxial pressure 
Pa = 4.5kbar along the o-axis.^^ Using those values, we find 7 = —31.6 degrees, r] = 0.92, A = 0.40, 
a = 1.18, and 0o = 32.6 degrees. Figure 4 shows the interlayer magnetoresistance for this parameter 
set. Prom the comparisons with the no tilt case (A = 1, a = 1.18) and the isotropic Fermi velocity case 
(A = 0.40, a = 1.0) with using the same other parameter values, we see that the ^-dependent interlayer 
magnetoresistance comes from the tilt of the Dirac cone. 

Experimentally a rough estimation of the tilt angle, 7, is obtained if we assume a ~ 1. In this case 
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Fig. 4. Dependence of the interlayer magnetoresistance on </> using the parameters estimated by Kobayashi et al. 
at Pa = 4.5kbar. Other parameters are set as -B = 6T, Bq = 0.5T, and 6 = 20 degrees. The results for no tilting 
case and isotropic Fermi velocity case are shown as well. 

anisotropy of the interlayer magnetoresistnace in the plane is mainly associated with the tilt. For a = 1, 
the function / takes a simple form, 



From this expression one can see that the resistivity takes the maximum in the direction of the tilt. An 
approximate value of 7 is found from the angle of that direction. The other parameters A and a can be 
estimated as follows. We first extract the function / {(p, a, 7, A) part from the experimental data using the 
parameters A and Bq determined from 6 dependence of p^^z ■ Taking the approximate values of 7 and A, 
which is obtained by using eq.(14), and a = 1 as initial values, more precise values are determined by the 
least squares method. 

Now wc comment on the conditions for the application of the formula (12). The formula is derived 
by using the zero-energy Landau level wave function. To justify this approximation, the magnetic field 
Bz should be large enough. The energy gap to the first excited Landau level is ei ~ 40\/S^K where 
Bz is measured in units of tesla if we assume the average Fermi velocity ~ lO^cm/s.^'^ Therefore, if the 
temperature is sufficiently lower than ei then the effect of the other Landau levels is safely neglected. As 
for the angle 9, the condition is 40%/ B sin > ksT. The formula including the effect of the other Landau 
levels is necessary for small Bz values. It is straightforward to extend the formula to this case. But the 
formula is complicated and such a formula is necessary only when one tries to see the crossover from the 
quantum limit to the semi-classical regime. However, this is beyond the scope of this paper. Another effect 
to be concerned is the Zeeman splitting. -"^^ However, the Zeeman splitting leads to a constant shift of the 
energy even though the shift depends on the spin. At fixed total magnetic field B the shift is unimportant 
for the determination of the parameters of the Dirac cone because it just leads to a modification of A in 



To conclude, we have derived the formula for the interlayer magnetoresistance in the presence of the 
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eq.(12). 
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tilt and the Fermi velocity anisotropy of the Dirae cone. The direction of the tilt is determined from 
the azimuthal angle dependence of the interlayer magnetoresistance. If the interlayer resistivity takes the 
maximum in some direction, then the Dirac cone is tilted in that direction. Physically the resistivity 
takes the maximum in the tilt direction because the interlayer current operator matrix element takes the 
minimum. The derived formula can be used to extract pressure dependence of the parameters of the tilted 
and anisotropic Dirac cone. It would be interesting to see the difference between the parameters determined 
by applying the formula to analyze the experimental data and the band calculation result. We expect that 
discrepancy between them should be large under low pressures because the Dirac fermions become unstable 
due to the electron correlation leading to the charge ordering^*^"^^ and/or superconductivity.^'^^ 
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